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Abstract. Relativistic particle production often requires the use of Tsallis statistics 
to account for the apparently power-like behavior of transverse momenta observed in 
the data even at a few GeV/c. In such an approach this behavior is attributed to 
some specific intrinsic fluctuations of the temperature T in the hadronizing system 
and is fully accounted by the nonextensivity parameter q. On the other hand, it was 
recently shown that similar power-law spectra can also be obtained by introducing some 
specific volume fluctuations, apparently without invoking the introduction of Tsallis 
statistics. We demonstrate that, in fact, when the total energy is kept constant, these 
volume fluctuations are equivalent to temperature fluctuations and can be derived 
from them. In addition, we show that fluctuations leading to multiparticle power-law 
Tsallis distributions introduce specific correlations between the considered particles. 
We then propose a possible way to distinguish the fluctuations in each event from 
those occurring from event-to-event. This could have applications in the analysis of 
high density events at LHC (and especially in ALICE). 
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1. Introduction 

Statistical modelling represents a widely used standard tool to analyze multiparticle 
production processes [1]. However, this approach does not account for the possible 
intrinsic nonstatistical fluctuations in the hadronizing system. These usually result in 
a characteristic power-like behavior of the single particle spectra and in the broadening 
of the corresponding multiplicity distributions. Such fluctuations are important as they 
can signal a possible phase transition(s) taking place in an hadronizing system [2], so to 
include such features, one should base this modelling on the Tsallis statistics [31 SJ [5] , 
leading to a Tsallis distribution, hq{E), which accounts for such situations by introducing 
in addition to the temperature T a new parameter, q > 1. This parameter is, as shown 
in [6l [7] , directly connected to fluctuations of the temperature. For g — )■ 1 one recovers 
the usual Boltzmann-Gibbs distribution, f{E): 



hq{E) = expg 





2-g 




T 






= ^ exp 1 





(1) 



hiE) = ^e.p^-^], (2) 

The most recent applications of this approach come from the PHENIX Collaboration 
at RHIC [8] and from CMS Collaboration at LHC [9] (see also a recent compilation 
[To]). One must admit at this point, that this approach is subjected to a rather hot 
debate of whether it is consistent with equilibrium thermodynamics or else it is only a 
handy way to a phenomenological description of some intrinsic fluctuations in the system 
^llj . However, as was recently demonstrated on general grounds in [12], fluctuation 
phenomena can be incorporated into a traditional presentation of thermodynamic and 
the Tsallis distribution (|2]) belongs to the class of general admissible distributions which 
satisfy thermodynamic consistency conditions. They are therefore a natural extension 
of the usual Boltzman-Gibbs canonical distribution i^!^. 

In fact, as was shown in |6j, assuming some simple diffusion picture as responsible for 
temperature equalization in a nonhomogeneous heat bath (in which local temperature, 
T, fluctuates from point to point around some equilibrium value, Tq) one gets the 
evolution of T in the form of a Langevin stochastic equation and distribution of 1/T, 
g{l/T), as a solution of the corresponding Fokker-Planck equation. It turns out that 
g{l/T) has the form of a gamma distribution, 

1 Tn / 1 Tn\^ / 1 Tn 
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I The recent generalization of classical tlicrmodynamics to a nonextensive case presented in [T5] 
should be noticed in this context as well. It is worth mentioning that, in addition to applications 
mentioned above, the nonextensive approach has also been applied to hydrodynamical models 14 and 
to investigations of dense nuclear matter, cf., for example, }15j . 
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Convoluting exp(— i?/T) with such g{\/T) one immediately gets Tsalhs distribution ([T]) 
with a well physically defined parameter q: according to Eq. ([3]) it is entirely given by 
the temperature fiuctuation pattern, which in turn is fully described by the parameters 
entering this basic diffusion process (like, for example, conductance and specific heat of 
the hadronic matter consisting this nonhomogeneous heat bath, cf., [6] for details). This 
approach was recently generalized to account for the possibility of transferring energy 
from/to a heat bath with a new parameter characterizing the corresponding viscosity 
entering into the definition of q (it appears to be important for AA applications [H [I6] 
and for cosmic ray physics [1 7J ; we shall not discuss this issue here) 1§|. 

In the next Section, we shall discuss the correspondence between fluctuations of 
volume V proposed in [20] and the presented above fluctuations of temperature T (both 
result in power-like distributions). Section |3] is devoted to a discussion of some speciflc 
g-induced correlations and fluctuations. Section H] is a summary. 

2. Fluctuations of V or T? 

We start by stressing that the form of g{l/T) as given by Eq. @ is not assumed, but 
derived from the properties of the underlying physical process in the nonhomogeneous 
heat bath. Apparently, the same results in what concerns the power-like character 
of single particle spectra and the broadening of the corresponding multiplicity 
distributions, P{N), were obtained in [20] without resorting to Tsallis statistics. It 
was assumed there that the volume V fluctuates in scale invariant way following the 
KNO form of P{N) deduced from the experiment |2l|j j1. We shall now demonstrate that, 
for the case of constant energy E considered in [20], both approaches are equivalent in 
the sense that one can start from fluctuations of V and recover fluctuations of T as 
discussed above, or else, one can start from fluctuations of T as given by g{l/T) and 
recover fluctuations of V as assumed in [20] . 

Following the approach of [20], for constant total energy, E = const, when both 
the volume V and temperature T are related via E ~ VT^, i.e., when 



§ In [T^ a similar suggestion of the extension of standard concept of statistical ensembles was 
proposed independently. A class of ensembles with extensive (rather than intensive) quantities 
fluctuating according to an externally given distribution was discussed. There is also a purely 
phenomenological approach treating the occurrence of a Tsallis distribution as a manifestation of the 
so called superstatistics, see [l9] . 

II It should be noticed that UA5 data [22] demonstrated that KNO scaling is broken via the energy 
dependence of the parameter k. In fact, as shown in [23], — —0.104 + 0.058 In . Therefore, in 
the scenario with fluctuations of the volume V, the scaling KNO form of the P{N) used to model these 
fluctuations is a rather rough simplification. On the contrary, in the scenario of the temperature T 
fluctuations, P{N) is given by a Negative Binomial Distribution, which adequately describes the data. 




(5) 
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the mean multiphcity in the microcanonical ensemble (MCE), N, can be written as 



This relation points to the KNO scaling form of the multiplicity distribution as good 
candidate for distribution of y, which was therefore assumed to be given by [2T] 

Hy) = Y^^y'''^M-ky). (7) 

From Eqs. and (E]) one gets that 



This means that T fluctuates as well according to the distribution ^/'((T)/T). The power- 
like form of the single particle spectra then follows immediately, all apparently without 
invoking any reference to Tsallis statistics. This completes the proposed picture which 
now comprises both fluctuations presented in the multiplicity distribution P{N) (from 
the scaling form of which one deduces the shape of volume fluctuations) and the power- 
like behavior of single particle spectra emerging because of temperature fluctuations that 
follow. Notice now that that ip{(r)/T) assumed here is, in fact, the same distribution 
as g{l/T) derived in Eq. (jl]) (with k = l/(g — 1), see also Eq. ( fT3ll below). So we obtain 
from V fluctuations the T fluctuations with the same functional form but now without 
the physical background behind Eq. (jl]) mentioned above. 

However, we can proceed in reverse order and obtain from T fluctuations (j4]) 
introduced in Section [T]the fluctuations of V introduced in [20], including the broadening 
of the corresponding P{N) which takes the form of a NB distribution. This point has 
been already shown in [21] and we shall quote here its main points for the sake of 
completeness. 

As was proved there, T fluctuations in the form of Eq. (jl]), discussed in Section 
m result in a specific broadening of the corresponding multiplicity distributions, P{N), 
which evolve from the poissonian form characteristic for exponential distributions to 
the negative binomial (NB) form observed for Tsallis distributions. One starts from the 
known fact that whenever we have independently produced secondaries with energies 
{i?j=i ...^at} taken from the exponential distribution f{E), cf. Eq. (j2]), i.e., when the 
corresponding joint distribution is given by 




where 






(9) 



and whenever 



N 



N+l 




(10) 



i=0 



i=0 



the corresponding multiplicity distribution is poissonian. 



P(iV) = ^ I exp (-A^) where N = —. 




(11) 
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On the other hand, whenever in some process particles with energies ... ^v} are 

distributed according to the joint A^-particle Tsalhs distribution, 

" - + l~N 

h ({-E'j=i,...,Ar}) = CAr 



A 



l-q 

(12) 



(for which the corresponding one particle Tsallis distribution function in Eq. ([T]) 
is marginal distribution), then, under the same condition ( ITU]) , the corresponding 
multiphcity distribution is the NB distribution, 

P{N) = J}^''\'ll, where k = (13) 

^ ^ r(Ar + l)r(A;) (jv>y^+'=)' q ' ^ ^ 



k J 

Notice that, in the limiting case of g — )■ 1, one has A; — ?■ oo and f|T3|) becomes a poissonian 
distribution (|TT]) . whereas for g — > 2 on has k 1 and f lT3|) becomes a geometrical 
distribution. It is easy to show that for large values of N and (A^) one obtains from Eq. 
( IT3|) its scaling form, 

/ N \ k'^ 

{N)P{N) = ^[' = J^^)= YXk)""' exp(-fcz), (14) 

in which one recognizes a particular expression of Koba-Nielsen-Olesen (KNO) scaling 
[2T] assumed in |20] to also describe the volume fluctuations, cf. Eq. This result 
closes the demonstration that, under the condition of constancy of total energy used 
here, T and V fluctuations are equivalent 0. 

We close this Section with the following remarks. As was said above, T fluctuations 
are derived from the more realistic description of the nonhomogeneous heat bath, and 
therefore parameter q and the Tsallis distribution, Eq. ([T]), reflect the physics of this 
heat bath. However, one can argue that, on this deeper level, this physics is nothing 
more than some phenomenological modelling, assumptions of which are reflected in q. 
The V fluctuations approach uses instead as its input the experimental knowledge of the 
scahng properties of particle multiphcity distributions, P{N), assuming that fluctuations 
presented there are transmitted to fluctuations of the volume. In fact, fluctuations of V 
could have some deeper phenomenological foundation, not mentioned in [20]. Namely, it 
is known that one observes experimentally a variation of the emitting radius (evaluated 
from the Bose-Einstein correlation analysis) with the charged multiplicity of the event, 
see, for example, [26]. An increase of about 10 % of the radius when the multiplicity 
increases from 10 to 40 charged hadrons in the final state was reported. Unfortunately, 

% In fact one can argue that the scahng form of P{N) visible in experiments points to the necessity 
of describing multiparticle production processes by means of Tsahis statistics. It is worth mentioning 
at this point that the connection between q and k was first discovered in [25| when fitting pp data for 
different energies by means of the Tsalhs formula ([1} . The resulting energy dependence of parameter q 
turned out to coincide with that of 1/fc of the respective NBD fits to corresponding P{N). It was then 
realized that fluctuations of N in the poissonian distribution p!T|) taken in the form of tp{N/ < N >), 
Eq. (HH), lead to the NB distribution 
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the quality of data does not allow us to precisely determine the power index of the 
volume dependence. It is also remarkable that both the energy density, pe = E /V , 
and particle density, pn = N/V, decrease for large multiplicity events. For N/ {N) ~ y 
one observes Pe/{pe) ~ and Pn/{pn) ~ y~^- AH these features deserve further 
consideration and should be checked in LHC experiments, especially in ALICE, which 
is dedicated to heavy ion collision. 



3. Some consequences of q statistics 

We would like to close with a short discussion of some consequences of g-statistics which 
can a priori be subjected to experimental verification: the g-induced correlations and 
event-by-event fluctuations. 

The g-induced correlations occur in a natural way in an A^-particle Tsallis 
distribution introduced in Eq. ( lT2l) . For Boltzmann-Gibbs statistics, for N 
independently produced particles, the joint distribution ([9]) can be written in factorizable 
form as a simple product of single particle distributions, 

N 
i=l 

However, such a product of single particle Tsallis distributions does not result in an A^- 
particle Tsallis distribution pi]- To get Eq. ( !T2|) one has to fluctuate the temperature 
in the distribution / ({-Ei=i,...,Ar}) above, i.e., 

POO 

h i{E,=^_^})= / / m=^_^}) g{l/T)d{l/T). (16) 

This procedure introduces correlations between particles. The corresponding covariance, 
Cov {Ei, Ej), and correlation coefficient, p, for energies are equal to 

Co.(E..fi.)= (3_^g-/_l3^) . (17) 
,^^-1. (18) 
As an illustrative example, we calculate the two-particle correlation function 

^'^^^"^'^-imm-y '''' 

This is shown in Fig. [T]for different values of variables 6 and a, 

Ei- Ej Ei + Ej 

6 = — and a = — (20) 



distribution, hg defined as 



Notwithstanding the rather complicated dependencies shown in Fig. [H the distance 
ig defined as 

POO 

hg{5)= / hgiE,,E, + 5)dE, (21) 
Jo 
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Figure 1. (Color online). Two-particle correlation function C2 {Ei, Ej) for Tsallis 
distribution with the parameter q = 1.125. 



for variables sampled from the joint distribution ([T^ is given by a Tsallis distribution 

W = ^[l-(l-?)5]^ (22) 

(analogously as the distance distribution for the exponentially distributed variables is 
given by exponential distribution, to which ( 12T|1 converges for g — > 1). 

In what concerns event-by-event fluctuations there are two types of fluctuations one 
can encounter, namely fluctuations from event to event or fluctuations in each event. 
Two scenarios are possible here. 

(i) In the first, T (and/or V) are constant in each event. However, because of different 
initial conditions, they fluctuate from event to event. In this case, in each event one 
should flnd an exponential dependence ([2]) with T = and a possible departure 
from it could occur only after averaging over all Nev events, k = 1,2,..., Nev It 
will reflect fluctuations originating in the different initial conditions collision from 
which a given event originates. Only inclusive distributions will be described by a 
Tsallis formula ([1]). Such a situation was discussed above. 

(ii) In the other scenario, T fluctuates in each event around some value Tq. In this case 
one should observe a departure from the exponential behavior already on the single 
event level and it should follow Eq. ([T]) with g > 1. This reflects the situation 
when, due to some intrinsically dynamical reasons, different parts of a given event 
can have different temperatures [H [3]. For volume fluctuations such a scenario 
seems to be unrealistic. 

In |27j we argued that event-by-event analysis of multiparticle production data are 
an ideal place to search for a possible fluctuation of the temperature characterizing 
a hadronizing source using the thermodynamical approach. Namely, an analysis of 
the transverse momentum spectra in Pb + Pb collisions at TeV energies of LHC 
(mainly in ALICE experiment) should allow us to distinguish between both scenarios 
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listed above. The LHC is designed for colliding proton-proton and nucleus-nucleus 
beams up to ^/s = 14 TeV |28]. Collisions at these unprecedented high energies will 
provide opportunities for new types of analysis. In proton-proton collisions at the 
highest possible energy, the expected charged particle multiplicity is only ~ 10 at the 
midrapidity region ( |?7| < 0.5) and it is roughly five times bigger in the full rapidity 
region |29] (the ALICE experiment has the possibility of measuring the distributions over 
the —5.0 < rj < 3.5 range, but the CMS and ATLAS experiments have a more limited 
coverage of |?7| < 2.5 units). The most important fact is that, in heavy ion collisions 
we have ~ A higher multiplicities (for central Pb + Pb collisions one expects ~ 2500 
particles at allowed rapidity acceptance region). This is enough to analyze event-by- 
event distributions over 3 orders of magnitude, which should allow us to determine the 
shape of the distribution in a single event. Moreover, in such circumstances we should 
be able to also construct the distribution dN/ d6 for 6 = pxi — Pxj for N{N — 1) = 6 ■ 10^ 
pairs in an event, over 6 orders of magnitude and using Eq. (12T|) test the above possible 
scenarios. Eq. (|2T]) tells us that, instead of distributions of pt, dN/dpr, one can use 
distributions dN/d6 and look whether it follows the Tsallis form on an event-by-event 
basis. Because for Eq. ( 12T|) one has A^(A^ — 1) entries to be compared with only for pt 
distributions, one expects that the distribution Eq. (I2T]) will reach further and it would 
be easier to differentiate between the Tsallis distribution and the usual exponential one. 




Figure 2. (Color online). The sensitivity to the parameter q of (a) the correlation 
function C2((T, (5) with a = 6 and (b) the distribution of 5, h{S), with error band 
corresponding to the event multiplicity N = 10'^. 

Fig. 12^ presents the sensitivity of the correlation function 02(0^,6) to different 
values of the parameter q chosen. Similarly, Fig. presents sensitivity of the proposed 
method in the event-by-event analysis on the example of distribution h{6) plotted for 
different values of the parameter q with error bands corresponding to the event with 
multiplicity A^ = 10^. 

We close this Section with a few remarks on what can be deduced from existing 
experimental data. To separate the two scenarios presented here one could, for example. 
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use the analysis of event-by-event fluctuations of the average transverse momenta pr 
presented at RHIC experiments. STAR data [30] show that Var {pj) / (pr)^ — 10"'', 
where pr denotes averaging in an event and {pr) is the overall event average. In the 
above estimation Var {px) = Var {j>T)data ~ ^'^^ (PT)mixed- STAR [30] measured also 
the quantity {ApTiApTj) = Cov {pTi,PTj) — Var (pr)- In this case ApTi = pn — (pt) 
was estimated without using mixed event and the resulting value of fluctuations was 
the same as above (cf. also [21] )■ Therefore, in the scenario of volume fluctuations 
(or temperature fluctuations occurring from event to event), assuming that pr oc T, 
we obtain Var{T)/{Tf = Var{pT) /{prf = 10-^ leading to g - 1 ^ 10"^ This 
is much smaller than the value g — 1 = 10^^ estimated from transverse momenta 
distributions [H |8l |9]. However, if the temperature T fluctuates within the event 
(for example, is different in every inter- nuclear collisions), we expect that for the Np 
projectile participants Var (pr) = Var{T) /Np. In this case, for Np = 100 the relative 
fluctuations are Var(T)/ (T)^ = 10"^, i.e., are comparable to those obtained from the 
transverse momentum distributions mentioned above. Similar estimations could be 
made also for [31] data. A more detailed analysis of this type is outside the scope of 
this paper and will be presented elsewhere. 

4. Summary 

To summarize: we have demonstrated that two approaches to fluctuation phenomena 
observed in multiparticle production processes as power law distributions or broadening 
of the corresponding multiplicity distributions, the one based on temperature 
fluctuations [6l [2l] and the one based on volume fluctuations [20] can be regarded, 
as long as the corresponding total energy is kept constant, to be equivalent. One 
can be deduced from the other. We have also shown that fluctuations which lead 
to multiparticle power-law distributions introduce some specific correlations between 
particles in the ensemble of particles considered and propose a way to distinguish the 
fluctuations in each event from those occurring from event-to-event by analyzing the 
distance distribution, Eq. (J2T1) . On the other hand, it seems that the already existing 
RHIC data on p^ correlations [301 EI] can also be useful to obtain such information (this 
point demands, however, a more detailed study). 
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